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We present more than 1000-day long photometry of EY Draconis in BV{RI)c passbands. The changes in the light 
curve are caused by the spottedness of the rotating surface. Modelling of the spotted surface shows that there are two 
large active regions present on the star on the opposite hemispheres. The evolution of the surface patterns suggests a flip- 
flop phenomenon. Using Fourier analysis, we detect a rotation period of Prot = 0.45875d, and an activity cycle with 
P ~ 350d, similar to the 11-year long cycle of the Sun. This cycle with its year-long period is the shortest one ever 
detected on active stars. Two bright flares are also detected and analysed. 
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1 Introduction 



Fast rotating late-type stars are known to possess strong 
magnetic fields (see e.g. lPizzolato et al 1[2003). These man- 
ifest in observable features like starspots, flares, emission in 
the Ha and Ca II lines, activity cycles, and significant emis- 
sion in_flie_^UV^X;rayregim e.g. Ayres & Li nskvl 
19801: iBerdyuginaL 120051: iHal 120081 IPagano et all I26oi1k 
Although Doppler-imaging and polarimetry gives a price- 
less tool in the hands of astronomers, the value of long-term 
photometric observations should not be underestimated, 
since this is the only way to continuously follow the be- 
haviour of a star on the timescale of years or decades. 



EY Dra is a well-known example of active stars with a 
history reaching back to almost two decades. In 1991 the 
ROSA T EUV/X-ray satelUte detected EY Dra as an EUV 
source. iJeffries et al. ( 1994 ) classified the object as a sin- 
gle, rapidly rotating (v sin i — 61kms^^) dMl-2e star, sim- 
ilar to HK Aqr. The authors presented a thorough spectro- 
scopic analysis and recommended the s tar as a candidate for 
studies of extreme magnetic activity. Eibd d 19981) carried 
out high-resolution Ha measurements, and suggested the 
presence of plage-l ike regions and prominence clouds above 
the stellar surface. Barnes & Collier Cameron used 
Doppler-imaging to obtain a surface map. The Doppler-map 
showed spots on all latitudes, but no significant evidence 
was found for the existence of a polar cap. Although ac- 
tive regions on the Sun appear always on lower latitudes, 
high-latitude spots and polar caps are familiar features on 
rapidly rotating stars as a result of the dominance of the 



Coriohs force over buoyanc y force on the rising flux tubes 
dSchussler & Solankil[l992l) . 

The first extensive photometri c study of EY Dra was 
done by iRobb & Cardinall (119951) . who examined V^-band 
data. The light curve showed an unusual W-shape, a 
small flare, and a period of 0.459d. Recent study of 



Korhonen et al. (l2007h used both V and R photometric 
measurements and spectroscopic data in optical and infrared 
domain. The authors find two active regions, chromospheric 
plages and prominences. The plages seemed to be a ssoci- 
ated w ith the active regions, as seen on the Sun. In I Vidal 
( l2007h F-band photometry was presented covering almost 
500d. The data indicated slow evolution of the surface spots 
and also a possible longer cycle of about 300d. The length 
of the dataset, however, was too short for such a conclusive 
statement. Plavchan et al.l (12009 ) analyzed 24 /im and 70 
fim Multiband Imaging Photometer for Spitzer (MIPS) ob- 
servations of 70 dwarfs, including EY Dra. They found no 
IR excess, and concluded there is no disk around the star. 

In this paper we carry out the analysis of BV{RI)c 
photometry of EY Dra covering more than 1 100 days. 

2 Observations 

Observations were obtained using the 60cm telescope of 
the Konkoly Observatory at Svabhegy, Budapest equipped 
with a Wright Instruments 750 x 1100 CCD camera (FoV 
17' X 24'). Measurements were carried out on 1 16 nights be- 
tween 2005 August 5 and 2008 October 24 using BV{RI)c 
filters. Altogether more than 2500 data points were col- 
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Fig. 1 Finding chart for the EY Dra field. 
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Fig. 3 Fourier spectrum of the EY Dra V light curve (top) 
and the spectral window (bottom). The insert in the top plot 
shows the power spectrum zoomed in to Prot ■ 
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Fig. 2 From top to bottom: Light curve of EY Draconis 
in V passband. The brighter points over the light curve at 
JDs 2453622 and 2454374 are results of flares. Dashed line 
shows where the data was split for Fourier analysis, as de- 
scribed in Section im Results of the spot modelling, show- 
ing the spot longitudes, sizes (A and 7, respectively) and 
spot temperatures. 



lected in each passband. Data reduction was carried out us- 
ing standard IRAF Q packages. Differential aperture pho- 
tometry was done using DAOPHOT package. 

GSC 03904-00259 and GSC 03904-00645 were used as 
comparison and check star, respectively (see finding chart in 
Fig.[T]l. Comparison-check magnitudes are —0.304 ±0.018, 
-0.316±0.014, -0.347±0.020 and -0.348±0.022 for B, 
V, Re, and Ic filters, respectively. Unfortunately, no stan- 
dard magnitudes are available for the comparison star, so we 
used standardized differential magnitudes in the paper. 

The resulting V light curve is plotted in the top panel of 
Fig. 12] For the phased light curves we used the ephemeris 
E = 2453588.16582 + 0.4587d. 



3 Analysis 
3.1 Period search 

Fourier analysis of the data was performed using M UFRAN 
(MUl ti frequency Analysis) written by iKollath 
(Il990h . The Fourier spectrum and the spectral window for 



the V light curve are plotted on Fig. [3] The strongest sig- 
nal is caused by the rotational modulation at Pi — Prot = 
0.45875 ± 0.00003d. The errors of the periods were esti- 
mated by increasing the residual scatter of the least-squares 
solutions by ™0005 whic h corresponds to 1 0% of the pre- 
cision of the data used (cf. IPress et al. . 1992 ). Another sig- 
nal at the half of the rotation period (P2 = 0.22936 ± 
O.OOOOld) is also present, and the result is the same when 
we fit the data with the main period and its first harmonic 
together This result reflects two major active regions on the 
stellar surface in the opposite hemispheres, (see the 'W- 
shaped light curves on Fig.|6]and the derived spot longitudes 



' IRAF is distributed by tlie National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science 
Foundation. 



© 2006 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



www.an-joumal.org 



Astron. Nachr. / AN (2006) 



791 



0.03 

0.02 

0.01 

0.00 

0.03 

.g 0.02 

E 0.01 
< 

0.00 
0.03 



iMi; ; JLLLl 



H 1 1 1 1 H 



H 1 1 1 1 H 



Fig. 4 Fourier analysis of the original V light curve split 
to three parts. The second row shows the Fourier-spectrum 
after prewhitening with Piot, in the third row both Piot and 
the signal at Prot/2 is removed. The same rotation period 
within the errors is found in all three parts. A signal at the 
half period is also present in all three segments. This sug- 
gests two major active regions on the stellar surface all the 
time. 
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Fig. 5 Pe riodogram ma de by using the SLLK method de- 
scribed by IClarkd (120021) . The large plot is zoomed in to 
Plot — 0.45875d, the insert shows the full periodogram in 
the range of 0. 1-0. 6d. 



from the modelling in Fig.lH After removing these signals 
plus a long-term trend which is comparable to the length 
of the dataset, we find a period of P3 = 348 ± 19d. Small 
residuals are present near the main period and its half, pos- 
sibly originating from the modulation residuals, i.e., from 
the changes of the light curve shapes. 

To check for changes in the periods, we divided the data 
to three parts, each covering « 400 days, as shown in the 
top panel of Fig. |2l and performed Fourier analysis for each 
segment. The result is plotted in Fig. |4] The resulting fre- 
quencies differ only a few seconds, and the differences are 
beyond the reliability of the period determination (O.OOOld), 
so the period is stable during the observations. 
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Fig. 8 Typical V light curves and colour index curves of 
BY Draconis showing two states of the surface with one and 
two large active regions. Continuous line shows the fitted 
spot model for V — Ic colour index. Note the tiny colour 
index amplitudes. 

Table 1 Basic stellar parameters for FY Dra 



spectral type 




dMl-2e ^ 


Tstar 




3900K2 


i 




f» 70°^ 


vsini 




61 km/s^ 


distance 




45.5 ±2.1 pc^ 


Mv 




8.54 ±0.12^ 


r sin i 




0.549 ± 0.002 Rq^ 




B 


V Rc Ic 


unspotted A mag. 


-0.470 


-1.325 -1.930 -2.600 


limb darkening^ 


0.763 


0.670 0.617 0.509 



Jeffries etal.1h994h 



Barnes & Collier Cameron! ( l200lh 
Hamme.(.1993E 



As a test, we checked the data for periods using the 
SLLK method (String/Rope l ength method using Lafler- 
Kinman statistic) described bv IClarke ( 2002 ). This method 
phases the light curves with different periods, and selects 
the period giving the smoothest light curve as the correct 
one. This method is very useful for finding periods for 
stars having non-sinusoidal light curve shape. The SLLK 
method gave the same result as the Fourier analysis (Prot = 
0.45875d, see the periodogram in Fig.|5]l. 



3.2 Spot modelling 

The light variation of FY Dra is caused by the rotation of the 
spotted surface. The modelling of this spo ttedness was done 
using SpotModeL dRibarik et all l2003b . which describes 
the intensity changes in an analytic way using homogeneous 
circular spots (Budding, 1977). This assumption is a valid 
approximation e.g. to the shape of the sunspots or the stellar 
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Fig. 6 Phased light curves of EY Draconis in V passband. These Ught curves served as input for the spot modelhng. 
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Fig. 7 Phased V — Ic colour index curves of EY Draconis. These curves were used to model the spot temperatures. 
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spots observed with Doppler-imaging. Another advantage 
of this model is the low number of free parameters: only 
the radius, location, and the temperature has to be fitted. 
"Spots" however are probably not homogeneous dark areas, 
but mixtures of unresolved dark (cool) and bright (hot) re- 
gions, similar to the active nests observed on the Sun. Thus, 
the temperature of the model spots is regarded as the aver- 
age temperature of hot and cool areas. 

During the modelling we have assumed two spots on 
the northern hemisphere, since a two-spot model can follow 
arbitrary spotted light curves wel l withi n the precision of 
the photometry dKovari & Bartusl Il997l) . SpotModeL is 
capable of spot temperature modelling making use of the 
colour index curves. We have chosen V — Ic (Fig. |7]l for 
determining the spot temperatures, since this colour index 
variation relates to the temperature changes the best. Dur- 
ing the fitting process T^a- was se t to 3900K according to 
Barnes & Colfier Cameron (l200lh Basic stellar parameters 
are summarized in Table[Tl The brightest values of the long- 
term light curve were used as unspotted magnitudes in each 
passband. Modelling was run for 16 light curves, plotted in 
Fig. |6] Two samples of the phased V light curve, B — V, 
V — Rc, and V ~ Ic colour indices are plotted in Fig. [8] 
together with the spot model fit. The fits show model curves 
with parameters derived from the modelling of the V — Ic 
colour index curve, only the temperature-dependent param- 
eters were changed for the plots. The resulting parameters 
are plotted in the lower panels of Fig. |2l except spot lat- 
itudes, since these values cannot be precisely determined 
at the given photometric precision. The spot temperature is 
about 500K below the photospheric temperature. Due to the 
low amplitude of the light curves the derived temperatures 
have relatively large errors, therefore their change seen in 
Fig.|2]has low significance. Models with different spot tem- 
peratures by 2-300K result in changing the spot radii by a 
few (1-3) degrees, so the overall picture of the spot radii 
variability remains the same. 

The longitudes of both of the two active regions remain 
similar, i.e. stay within « 80° during the time of the obser- 
vations, more than 1000 days. In this time, the active region 
at phase « 0.3 slowly decays, while the other one, at phase 
« 0.8 gets more and more prominent. 

3.3 Flares 

During the observations two noteworthy flare events hap- 
pened: in JDs 2453622 and 2454374, lasting about 3.5-3.8 
hours. Both events could be well seen in all filters, as plotted 
in Fig. |9l The plots show instrumental light curves, which 
was used for flare energy calculations. The fast colour in- 
dex changes does not allow a perfect transformation to the 
international system during flares; the necessary interpola- 
tions (which would be only rough approximations of the 
colour indices anyway) would blur the details of the ob- 
servations. The first event is a complex one, covering three 
days. Before the large eruption five smaller outbursts were 
observed at different phases, probably not connected with 



Table 2 Estimated energies, equivalent durations of the 
two large flares, and quiescent fluxes in different filters. 





Flare 


Equivalent 


Quiescent 




energy 


duration 


flux 




[xlO^'' 


ergs] 


[sec 




[xlO^^ ergs/sec] 


Flare 


#1 


#2 


#1 


#2 




B 


4.53 


4.82 


1232 


1311 


3.68 


V 


1.77 


2.13 


292 


352 


6.06 


Rc 


2.97 


3.35 


206 


232 


144 


Ic 


1.01 


2.70 


82 


219 


12.3 



#1: JD2453622, #2: JD2454374 

Table 3 Values for flare temperature Tp and relative area 
A = Ap/Astar as a result of our trial-and-error procedure 
(see text). 



Flare 


#1 




#2 




Tflaro 


Ab 


Av 


Ab 


Av 


[K] 




[ 






8000 


1.04 


0.29 


0.56 


0.13 


12000 


0.25 


0.09 


0.13 


0.04 


16000 


0.12 


0.04 


0.06 


0.02 


20000 


0.07 


0.03 


0.04 


0.01 


25000 


0.05 


0.02 


0.03 


0.01 



#1: JD2453622, #2: JD2454374 

the main event. All of them, however happened when the 
dominant active region was visible, around mid-phase of 
0.3, suggesting a connection between the starspot and the 
flares. According to its phase of about 0.9, (cf. Fig. |6]and 
|2]l, the other flare on JD 2454374 could also be associated 
with the smaller active region on BY Dra. 

For estim ating flare energies we followed the method 
described by iKovari et alj (120071) . The energy emitted by 
the two flares and their equivalent durations are summarized 
in Tab. |2l The second flare (JD 2454374) is slightly more 
energetic, though its peaks are lower in the red colours than 
that of the first flare (JD 2453622) (Fig.|9t>). 

We give a rough estimation for the colour temperature of 
the two flares at thei r maxima, by following the method of 
de Jager et al. Il986h . The flare flux is significantly measur- 



able in B, and barely in V, thus we use only those colours. 
If A is the ratio between the projected area of the flaring 
region and the visible stellar surface, then the flux ratio be- 
tween the flaring and the quiescent star is 



[Fs(l -A) + FfA]/Fs = 1-A + ABf/Bs 



(1) 



where Fg is the total flux of the quiescent star (Table|2]i, Fp 
is the flare flux. Bp and Bs are the blackbody fluxes of the 
flare and the star, respectively. Note that there are two un- 
known quantities in Eq.[Tl A and Tp, the flare temperature 
through Bp. We only search for solutions by trial-and-error 
in a reasonable temperature range (see Table [3]). If we as- 
sume the given range feasible for a dMe star, we get some 
0.1-0.01% of the visible stellar surface for the flaring re- 



gion, which is comparable with other results (cf .Olah et alJ 
200lh . 
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Fig. 9 (a): Phased bv light curves showing the three nights of the first flare event separately. Vertical arrows show the 
outbursts, dashed lines and the horizontal arrows at the bottom of the plots indicate the beginning of the observations, (b): 
Phased bv{ri)c light curves showing two energetic flares on EY Dra. Gray points show additional observations of nearby 
nights for comparison, v, rc, and ic light curves are shifted arbitrarily. Note, that the plots show instrumental values (see 
text). 
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Fig. 10 Histograms showing the longitude distribution of 
the two active regions (a), and the longitude difference be- 
tween the spot centers (b). 



4 Discussion 

Our spot modelling revealed two permanent, large active re- 
gions on EY Dra. From the beginning of the observations to 
JD «2454600 we can follow the decay of the spot group 
around phase 0.4 (see spot radii values on Fig.|2]and phased 
light curves on Fig. |6l). At the same time, the active region 
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Fig. 11 V light curve showing the w 350d period after 
prewhitening with the rotational modulation and the long- 
term trend. 



around phase 0.8 became more prominent. During the ob- 
servations the total spotted area decreases from w 15% to 
« 10%. Note, that the derived spottedness values are lower 
limits, since the unspotted brightness is set to the maximum 
observed magnitudes. The change of the spot parameters 
in Fig, m i s remin iscent to those observed on FKCom by 
Olah et alJ (l2006l) . The longitudes of the spots remain in a 
relatively narrow region, the distance of the two active re- 
gions is also similar, « 130° (see the two histograms in 
Fig. [TOb . The third panel of Fig. |2] shows clearly, that af- 
ter JD 2454400 the dominance between the two active re- 
gion exchanges; the light curves of Fig. |6] directly show 
this change (JD 54363-54377 and JD 54539-54567). This 
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kind of so-called flip-flop phenomenon was observed on 
many stars, see e.g. lJetsu et alj (ll99ll) . lElstner & Korhonen 
(l2005h showed that a mixed-mode dynamo including a non- 
axisymmetric and an oscillating axisymmetric mode can de- 
scribe this kind of phenomenon in stars with thick convec- 
tion zones. 

The Fourier analysis showed, that there is one strong 
signal in the light curve. This modulation with Prot ~ 
0.45875d is caused by the rotation of the spotted surface. 
Since in most cases there are two active regions on EY Dra 
(see Fig.|6]l, a signal appears also at Piot/2 in the spectrum. 
Next to the strong signals of P^ot and Prot/2 weaker peaks 
can be found. This might be the result of differential ro- 
tation, but the difference between the signals is too small 
to draw such a conclusion. Differential rotation can still be 
present on the stellar surface, but when the latitudes of the 
active regions are close to each other, it is beyond the limit 
of detection. Unfortunately from photometry alone no trust- 
worthy information can be determined on the spot latitudes. 
On the other hand, the small peaks near the rotational period 
and its half in the Fourier spectrum may originate simply 
from the modulation residuals, i.e., from the change of the 
light curve shapes. In summary: the detected rotational pe- 
riod of EY Dra is remarkably stable during the three years of 
our observations. This could be the signature of a persisting 
magnetic configuration on the stellar surface, with spots on 
the same latitude. In the same time the two narrow regions 
of spot longitudes (see Fig.fTOh) also means long-term sta- 
bility of the spot position. Stable magnetic fiel ds lasting for 



at lea st one year was observed on V374 Peg dMorin et al. 



20081), which is about 0.3 solar mass and fully convective. 
EY Dra is also a low mass star (about 0.5 solar mass, cf. 
EibeL Il998h with a very deep convection zone, therefore a 
similar stable magnetic configuration on its surface is pos- 
sible. 

Another peak in the amplitude spectrum indicates a 
long-term variation with Pcycie ~ 350d (see the light curve 
plotted in Fig.fTTI). This can be a result of an activity cycle, 
similar to the 1 1-year long cycle observed on the Sun. Signs 
for the presence of this year-long c ycle can als o be found in 
ROTSE archive data, as shown in I Vidal (120071) . This solar- 
like cycle is the shortest activity cycle known. As to our 
knowledge, there are only two M stars which show activity 
cycles: EY Dra (Ml-2) with a very fast rotation of 0.46 days 
and a cycle length of about 1 year, and HD 95735 (M2), 
slowly rotating with a period of about 55 days hav ing cy- 
cles of 3.4-3.9 and «11 years (cf. lOlah et aP (I2006h . 



5 Summary 

- Fourier analysis showed, that the rotation period is 
-Prot — 0.45875d, which is stable during the observa- 
tions. 

- An activity cycle of « 35Gd has been found. 

- Spot modelling showed that the active regions on the 
surface are located around « 130 and « 250 degrees. 



- The evolution of the surface indicates a possible flip-flop 
mechanism. 

- No reliable evidence was observed on differential rota- 
tion. 

- Two flare events were found around JDs 2453622 and 
on 2454374, possibly connected with the active regions 
on the surface. 
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